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ABSTRACT 

Submillimeter cameras now have up to 10 4 pixels (SCUBA 2). The proposed CCAT 25-meter submillimeter 
telescope will feature a 1 degree field-of-view. Populating the focal plane at 350 microns would require more than 
10 6 photon -noise limited pixels. To ultimately achieve this scaling, simple detectors and high-density multiplexing 
are essential. We are addressing this long-term challenge through the development of frequency-multiplexed 
superconducting microresonator detector arrays. These arrays use lumped-element, direct-absorption resonators 
patterned from titanium nitride films. We will discuss our progress toward constructing a scalable 350 micron 
pathfinder instrument focusing on fabrication simplicity, multiplexing density, and ultimately a low per-pixel 
cost. 
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1. INTRODUCTION 

Modern cryogenic detector arrays for infrared and submillimeter astronomy are now routinely achieving the 
kilopixel scaleP-"-^ While this represents a significant technological milestone, there remains a strong demand for 
even larger arrays. For example, the proposed CCAT 25-meter submillimeter telescope^ would require more than 
10 6 photon-noise limited pixels in order to fully sample the focal plane at a wavelength of A = 350 /.tm. In order to 
achieve this considerable scaling, candidate cryogenic detector technologies must be highly multiplexed, simple 
to fabricate, and the per-pixel cost must be minimized. One particularly attractive technology for achieving 
very large arrays is the superconducting microresonator detectorPtH The pixel in this case consists of a radio- 
frequency resonator patterned from a superconducting thin film. Incident energy absorbed by the resonator 
breaks superconducting Cooper pairs generating quasiparticles. This perturbs the complex conductivity of the 
thin film 6a(u>) = 8cr\ — jScri, consequently producing a measurable shift in the frequency io r and quality factor 
Q r of the resonance. Strong frequency-domain multiplexing is achieved by coupling an entire array of pixels to a 
radio-frequency transmission circuit; the frequency of each pixel is tuned by adjusting the resonator reactance. 

Two methods are typically employed to couple optical signals into superconducting microresonator pixels. 
One approach is to collect the incident radiation with an antenna and deposit the energy in the high-current 
portion of a quarter-wave coplanar waveguide resonator PQ21 This method has a number of advantages including 
good beam control, frequency selectivity, and the possibility of multi-chroic operation through the use of on-chip 
bandpass filters. The second method utilizes a lumped-element resonator consisting of an inductive meander 
section with an attached intcrdigitatcd capacitor and is commonly referred to as the lumped-element kinetic 
inductance detector (LEKlD)! 13 ^ 1 ' In this design, radiation is directly absorbed by the inductive meander which 
is carefully designed to provide an effective far-infrared surface resistance of around i? c ff = 377 f2/(l + ^/e7) ~ 86 
when back-side illuminating the pixel through a silicon substrate with dielectric constant e r 11.5. This result 
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Figure 1. The MAKO instrument, (a) 3 He measurement cryostat with a stable base temperature of 250 mK. (b) First 
generation 100 pixel test array. Here the single transimssion line is coupled to the inductive meander portion of each pixel 
and is used to read out the entire array. For this initial design the capacitor occupies ~ 35% of the pixel area. 

is modified to R c s = 377 £1/^/1^ if a backshort is used. An important advantage of the lumped-element approach 
is that radiation absorbed by the meander is uniformly distributed over the inductor volume. This avoids regions 
of high quasiparticle density which reduces the quasiparticle lifetime and consequently the device responsivity. 
Lumped-element designs are also typically less complex to fabricate than antenna-coupled pixels, often requiring 
at most a few layers of lithography to realize an entire array. On the other hand, good baffling and stray light 
rejection is relatively more important for absorber-coupled detectors. 

For the LEKID, achieving a far-infrared surface resistance of i? c g s» 86 f2 constrains the sheet resistance R s 
and area filling factor t\a of the superconducting film by i? c ff ~ Rs/i]A- To present a uniform absorber to incident 
radiation of wavelength A, the meander spacing is further constrained by s -C A. For A = 2 mm, near background- 
limited performance has been demonstrated using aluminum LEKIDs for ground-based astronomy.^ In this case, 
a low filling fraction T)a was used to compensate for the low resistivity of aluminum. At shorter wavelengths, 
the required meander spacing precludes the use of traditional elemental superconductors for reasonable film 
thicknesses and linewidths. Fortunately, the discovery by Leduc et alP^ of the exceptional surface quality Q a of 
superconducting titanium nitride (TiN) provides an alternative. In this case, the resistivity pxtN ~ 100 fj.fl/ cm 
is several orders of magnitude higher than standard metals, allowing high fill fractions to be used with reasonable 
film thickness and optically lithographed linewidths. Complete LEKID design details for the submillimeter/far- 
infrared are presented in an article in these proceedings by McKenney et alP^ 

We are now actively engaged in constructing MAKO, a scalable 350 micron pathfinder instrument emphasizing 
multiplexing density, fabrication simplicity, and a low per-pixel cost. An image of the cryostat and a 100 pixel 
first-generation test array are shown in Fig. [3l We have chosen LEKIDs patterned from TiN films for this project 
as the most suitable detector for achieving the project goals in the far-infrared. Notably, these pixels are designed 
to resonate at f=a 100 MHz frequencies. All previous kinetic-inductance based imaging arrays have operated at 
> 1 GHz. This shift toward lower frequencies provides numerous advantages including an improved pixel noise 
equivalent power (NEP), a simplified analog readout circuit, and a higher achievable multiplexing density. These 
low frequency advantages are discussed below along with a detailed description of the cryo-mechanical and 
readout systems currently being commissioned. 



2. LOW FREQUENCY ADVANTAGE 



The NEP of a superconducting microresonator detector is given by: 
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where Eq. ([1} and Eq. @ are the NEP when measuring perturbations of the resonator dissipation (Qr 1 ) and 
frequency, respectively. In both of these equations the first three terms are the photon, generation-recombination, 
and amplifier noise contributions. The full expression for the generation-recombination noise can be found in a 
recent review by Zmuidzinas 8 ^ but is not a significant noise source for the current designs where hv 3> A. 

The additional term in Eq. ([2]) reflects frequency noise due to dielectric fluctuations in the resonator capacitor. 
Commonly referred to as two-level system (TLS) noise, these fluctuations arise from surface oxides or other 
amorphous dielectric material on the surface of the substrate! 17 ^ 18 ! TLS noise has been extensively studied and a 
variety of techniques have been developed to mitigate this effect CsHU] One of the simplest methods for reducing 
TLS noise while simultaneously reducing the amplifier noise contribution is to increase the readout power. For 
the simplest readout protocols, readout power is limited, this is limited by the onset of nonlinear effects and the 
eventual bifurcation of the device response. Bifurcation occurs when 



where Pg is the generator power at bifurcation, Q c is the resonator coupling quality factor, a c = .77 is the 
nonlinearity parameter at bifurcation, and E* sets the scaling of the nonlinearityP^ If the bifurcation mechanism 
is due to the nonlinearity of the kinetic inductance, E* is related to the condensation energy of the inductor 
E CO nd = NqVl& /2, where Nq is the material-dependent single-spin density of states at the Fermi energy, 
A fa 3.5fcsr c /2 is the superconducting gap, and Vl is the volume of the inductor For TiN, N ps 8.7 x 10 9 
/im~ 3 eV _1 P^ Realizing high multiplexing densities requires a high resonator quality factor Q r . According to 
Eq. (|3l), this lowers the onset of bifurcation and consequently degrades the NEP when operating below this 
limit. This can be avoided by decreasing Q r however this increases the resonator bandwidth and consequently 
the per-pixel cost. While it has been shown for a TiN far-infrared LEKID pixel that it is possible to operate 
at powers beyond the onset of bifurcation,^ this requires a considerable increase in readout complexity which 
inevitably leads to increased pixel cost. 

Eq. @ suggests another method for improving the device NEP while operating below the onset of bifurcation. 
The ratio of the reactive to dissipative response is given by (3 — Saz/Sax ps 2kT/hf. Historically, kinetic- 
inductance detectors have been designed to operate around a few GHz. Operation in this regime is useful for 
providing the large per-pixel bandwidth necessary for single-event detection such as single-photon counting or 
phonon sensing J^HUD However, the fundamental bandwidth requirement for submillimeter-wave astronomical 
imaging is set by the telescope scan speed. A per-pixel bandwidth of ~ 100-1000 Hz is sufficient for this 
application. If for multiplexing reasons we require Q r to be of order 10 5 , this sets a lower limit on the resonance 
frequency of f r ~ 10-100 MHz. Designing the pixels to resonate at MHz instead of GHz frequencies increases /3 
by more than an order of magnitude with a proportional improvement in NEPf req . 

Operating at lower frequency provides other crucial advantages apart from improving NEPy re(? . Readout 
for the GHz regime typically begins by generating a frequency comb of measurement tones with a spacing 
corresponding to the pixel frequencies using a high resolution digital-analog converter (DAC). These tones are 
located in a few hundred MHz of bandwidth around DC. The tones are then shifted to the desired GHz band 
using an up-converting mixer and a separate local oscillator. After interacting with the resonances the waveform 
is then shifted back down to baseband by a second mixing stage. Subsequent digitization is accomplished using a 
high resolution analog-digital converter (ADC). By operating below a few hundred MHz, the two analog mixing 
steps can be removed. This reduces readout cost and complexity as well as eliminating analog mixer distortion. 
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Figure 2. Collision probability per octave of readout bandwidth. In this calculation, a fixed number of pixels with 
resonator quality factor Q r are evenly distributed throughout an octave of bandwidth. Their positions are then perturbed 
randomly from their initial value /,• with a gaussian distribution of width Uf r / ,f r - Two pixels are assumed to collide (ie 
suffer from unrecoverable electronic cross-talk) if their operating frequencies are within five linewidths. In this case, both 
pixels are removed. Typical values from previous resonator detector arrays of Oj T j j T — 2 x 10 -3 (solid) and 5 x 10~ 4 
(dashed) are shown. Note that as the Q r increases, the linewidth and consequently the collision probability decrease for a 
fixed number of pixels per octave. If readout cost dominates, potentially it is reasonable to accept a higher percent pixel 
loss if the total number of pixels per readout increases. 

The available electronic readout bandwidth is usually constrained by the sampling speed of the high resolution 
(> 12 bit) ADC. Currently this provides at most a few hundred MHz of readout bandwidth. By moving from 
GHz to MHz frequencies, the per-pixel bandwidth A/ = f r /Q r decreases allowing a greater pixel count per fixed 
electronic bandwidth. A statistical calculation of the expected pixel loss due to resonance collision for various 
packing densities is shown in Fig. [2] Here, the fractional error in the resonator design frequency is assumed to 
be set by lithographic constraints. The series of curves correspond to packing various pixel counts per octave of 
bandwidth. As the resonator Q r increases, fewer pixel collisions occur as a result of the decreased per-resonator 
bandwidth. The important point here is that the number of octaves per fixed electronic bandwidth increases at 
lower frequencies. For example, 250 MHz of electronic bandwidth only occupies about a quarter octave in the 
range 1-2 GHz whereas the same electronic bandwidth spans more than two octaves from 50 to 250 MHz. As 
readout cost is roughly proportional to electronics bandwidth, increasing the multiplexing density directly lowers 
per-pixel cost. 

3. MEASUREMENT SYSTEM 

3.1 Optical Design 

For simplicity, MAKO has been design to use the same relay optics system as the existing SHARC-2 instrumenlP3 
at the CSO. The instrument will be mounted at one Nasmyth focus of the telescope, with two flat fold mirrors 
relaying the f/12.96 Cassegrain focus of the telescope to a single off-axis elliptical mirror reducing the focal ratio 
to f/4.48. The tertiary image of the far field, 45 mm in diameter, is formed 320 mm from the elliptical mirror 
with a 32 mm diameter image of the primary mirror at 178 mm. These positions lie 193 mm and 51 mm from 




Figure 3. Cryomechanical and optical layout of the MAKO instrument. This system is designed as a drop-in replacement 
for the SHARC-II instrument currently installed at the Caltech Submillimeter Observatory. 



the cryostat reference plane, respectively. The final focus offers a 3.7 arcminute field of view, limited by the tube 
through the telescope elevation bearing. 

The cryostat window is composed of high-density polyethylene, 1 mm thick with a clear aperture of 70 mm. 
Long-wavelength-pass filters are mounted at the apertures of the 50 K and 4 K radiation shields in order to 
reject infrared radiation from the cold stages of the instrument. At the 50 K shield, a Z-cut crystal quartz 
filter of diameter 65 mm and 2 mm thickness, with a 58 /xm clear low-density polyethylene (LDPE) coating 
applied to the upper surface and a 58 /im black LDPE layer adhered to the lower surface. The LDPE layers act 
as an anti-reflection coating for the quartz filter, while the black LDPE is also useful for blocking optical and 
near-infrared radiation. As second, identical quartz filter is mounted at the aperture of the 4 K radiation shield, 
in addition to a 100 /im and a 300 fim metal mesh long-wavelength-pass filter, produced by QMC Instruments^] 
The observing band of the detectors is defined by a further metal mesh bandpass filter, also manufactured by 
QMC Instruments, with a band center at 340 and a 10% bandwidth. The bandpass filter is integrated with 
the detector package. The peak transmission of the full filter system is approximately 50%. 

3.2 Cryomechanical Design 

Since MAKO is designed to be a "drop- in" replacement for SHARC-2, it was necessary that the new cryostat be 
mechanically compatible with the existing hardware. Unlike SHARC-2, the MAKO cryostat uses a Cryomech 
PT4-10 pulse tube mechanical cooler to provide the primary heat uplift, whereas the current SHARC-2 system 
is a liquid nitrogen/liquid helium system. The pulse tube provides a cooling capacity of 40 W at 45 K on the 
first stage, used to cool a radiation shield and optical filters as well as intercepting the conduction load from 
wiring, while providing a cooling capacity of 0.9W at 4K on the second stage. The 4 K stage of the cryostat 
cools a second radiation shield with further filters and provides the heat sink for the 3 He sorption refrigerator 
that cools the detector array (described below). The vacuum jacket and radiation shields, shown in Fig. [3] are 
an all-aluminium welded design, fabricated by Precision Cryogenics IncQ The outer shell of the cryostat has a 
diameter of 470 mm and is 1250 mm in height, including the open frame protecting the top plate. 

The 4 K space of the cryostat is 290 mm in diameter and 420 mm in height. In addition to the outer 4 K 
radiation shield, an inner shield, composed of Aluminium 6061, encloses the sub-Kelvin stages of the cryostat. 
This inner shield provides additional baffling against stray light and the mechanical mounting of a cold stop 
aperture at the image of the primary mirror. This structure is mounted to the 4 K cold plate via a CFRP 
hexapod structure. Although not intended to be thermally isolating, CFRP represents an excellent structural 
material being inexpensive, light and extremely stiff. Cooling is provided by a separate copper strap, with a 
reach-around strap cooling the cold stop aperture itself rather than relying on conduction through the aluminium 
structure. 

The detector package is cooled by a three-stage helium sorption fridge, design and fabricated by Chase 
Cryogenics @ with a 4 He cold head and two 3 He cold heads. Unlike other designs of sorption cooler in which the 
4 He stage is integrated with the first 3 He cold head, the MAKO fridge has three separate cold heads, although 
the 4 He head is used only as the condensation point for the 3 He in the other heads. The first 3 He cold head 
cools an intercept stage to approximately 350 mK, buffering the conduction load from the 4 K stage and the 
superconducting coaxial cables running to the detector package. The detector package is cooled to 240 mK by 
the second 3 He head. The intercept and detector stages are thermally isolated with CFRP bipod trusses and 
are connected to the sorption cooler with compliant copper foil straps. 

3.3 Electronics 

State-of-the-art readout electronics for microresonator detectors currently achieve multiplexing factors around a 
couple hundred, with a record of 256!^^ As these electronics currently cost several thousand dollars, the per- 
pixel cost is many tens to upwards of a hundred dollars. Simply scaling these existing systems, the readout cost 
will completely dominate any system achieving tens of kilopixels. Operating at lower frequencies results in narrow 
pixel linewidths which necessitates increased frequency resolution for both the individual tone generation and 

*http:/ /www. terahertz. co.uk 
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post-transmission signal separation compared with existing readouts. However as already discussed, the narrow 
linewidths provide more than an order of magnitude gain in multiplexing density which directly translates to 
decreased pixel cost. In order to realize this significant gain, we are meeting the challenge of high-resolution read 
out with the two electronics solutions discussed below. 

Our first readout places a high priority on minimizing development time. This approach uses a dual-channel 
commercial digitizing card (Pentek, model 78650). This card is equipped with a pair of DACs each connected 
to a circular buffer enabling looped playback of a precalculated waveform. Dual ADCs on the card are then 
used to digitize the transmitted waveforms. Rather than processing this received waveform locally, the data is 
streamed via a PCIe communication bus to a graphical processing unit (GPU). The GPU performs a continuous 
FFT on the received data. Selection of only the relevant FFT spectral bins and averaging then reduces the 
data rate before saving the results to disk. While a field-programmable gate array (FPGA) is included on the 
card, it only is used for communication and signal routing and does not perform any digital-signal processing. 
All FPGA functionality is preprogrammed and accessed using standard CH — h libraries. No specialized hardware 
development language (HDL) programming is required. 

The most challenging aspect of this approach is processing the data provided by the digitizing card in real 
time. The dual 12-bit ADCs each sample at a rate of 500 megasamples per second (MSPS). The PCIe bus 
provides 500 MB/s per lane of bandwidth and the digitizing card utilizes 8 lanes, thus the 4 GB/s of total 
throughput comfortably supports the 1.5 GB/s produced by the ADCs. Once the data has been received by 
the GPU, it must compute successive N point FFTs on the fly. For an ADC sampling at 500 MSPS, an FFT 
of at least 2 22 points must be continuously calculated to achieve a ~ 100 Hz resolution. The computational 
cost of an FFT is « 5N \og 2 (N) floating-point operations (FLOP) and thus the rate the gpu must support is: 
(ADC Sample Rate/iV)*5iVlog 2 (7V) FLOP/s. For N between 2 20 to 2 26 , this corresponds roughly to 100 to 
130 GFLOP/s to continuously support both ADCs at their full sample rate. Modern GPUs are benchmarked 
to provide more than 150 GFLOP/s for single-precision FFTs providing a reasonable margin for computational 
overhead. 

The second readout we are currently developing places a premium on algorithm efficiency and hardware 
cost reduction. Currently in the initial design phase, this readout will utilize an FPGA-based ROACH plat- 
form to evaluate potential signal-processing algorithms.''' 5 A variety of techniques, including multistage FFTs, a 
polyphase filter-bank, or a polyphase filter-bank followed by a second stage of digital down-conversion as demon- 
strated by McHugh et alP^l are being considered. A portion of the signal processing may be performed by a host 
computer after an initial stage of processing reduces the information rate sufficiently for communication transfer. 
The balance between FPGA and CPU/GPU processing will be optimized for cost and simplicity. Eventually, we 
plan on migrating this readout to custom hardware!^ Applying these optimizations should result in a readout 
cost reduction of around an order of magnitude. Ultimately, the electronics will likely feature an application 
specific integrated chip (ASIC) in order to reduce the volume cost and readout power requirements crucial for 
remote-location and space-based deployment of large format instruments. 

Comparison of the measurement signal to the amplifier noise affords an estimate of the highest multiplexing 
factor achievable on a single measurement line. As discussed above, the maximum achievable signal power is 
limited by the onset of bifurcation. For our current pixel design, Vl = 11550/im 3 , T c ~ 1.2K, Q c = 7 x 10 5 , and 
Q r ~ 5 x 10 4 under 100 pW (single-polarization) of illumination from a blackbody. Using Eq. ([3]) and setting 

(c) 

E* = Econd, the predicted bifurcation readout power is Pg w 307 fW = -95 dBm. Measurements of our actual 

(c) 

pixels however have found an onset of bifurcation at a significantly higher Pg — 3.2 pW = -85 dbm. This is the 
incident microwave power from the generator. Only a portion of this is transmitted to the amplifier; the rest is 
absorbed by the pixel or reflected. For optimal coupling (Qi = Q c ; % c = 1) a quarter of Pg C ^ is transmitted to 
the amplifier. Thus at bifurcation, the signal at the amplifier is of order 800 fW. A noise temperature of T a = 4 
K is typical of widely available cryogenic silicon-germanium bipolar transistor or HEMT amplifiers. The noise 
spectral density of the analog electronics is then calculated to be 101og 10 (/cT a /P tra „ s ) = —101 dBc. For both 
the readout options discussed above, the ADC is a Texas Instruments ADS5463. The SNR for this ADC is ~ 56 
dBFS. Assuming white noise and a sample rate of 500 MSPS, this corresponds to a noise spectral density of -143 
dBc/Hz. Comparing these results, around 10 4 tones can be simultaneously multiplexed before the ADC begins 
to degrade the measurement signal-to-noise ratio. 



4. CONCLUSION 



MAKO is a 350 micron camera currently being constructed and tested at Caltech and JPL. An on-sky demonstra- 
tion of MAKO at the Caltech Submillimcter Observatory is planned for late 2012. By emphasizing multiplexing 
density, fabrication simplicity, and a low per-pixel cost we aim to demonstrate that economical, large-format, 
far-infrared imaging arrays are readily achievable using superconducting resonator detectors. 
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